Detection of vorticity in Bose-Einstein condensed gases by matter-wave interference 
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A phase-slip in the fringes of an interference pattern is an 
unmistakable characteristic of vorticity. We show dramatic 
two-dimensional simulations of interference between expand- 
ing condensate clouds with and without vorticity. In this way, 
vortices may be detected even when the core itself cannot be 
resolved. 

Bose-Einstein condensation of a dilute gas has now 
been accomplished for three different alkali atoms by a 
number of research groups . Many of the properties of 
the gas which were calculated from mean-field theories, 
such as condensate density, collective excitation spectra, 
and the speed of sound have also been found to agree 
well with these theories. One of the untested predictions 
of mean-field theory is the quantization of circulation of 
the fluid, leading to quantized vortex states. In this Let- 
ter we show how such a state may be detected by using 
matter-wave interference. This method may allow us to 
follow the dynamics of states with vorticity as well. 

Several ways of creating vortices in a Bose-Einstein 
condensate (BEG) have been suggested. Initial experi- 
mental attempts relying on a rotating, off-resonance laser 
beam were inconclusive Another idea is that vor- 
tices may appear spontaneously out of fluctuations, as 
a result of rapid cooling through the critical tempera- 
ture @. Both of the fore going methods are analogous to 
techniques in more traditional superfluids such as liquid 
^iJe. Unlike those systems, the internal states of the al- 
kali atoms provide a convenient "handle" for the optical 
creation and detection of vortices. Starting with the trap 
ground state, vortex states can be reached via a Raman 
transition quickly or adiabatically ||]. 

Here we focus on the detection of states with vorticity 
or non-zero winding number. Indirect detection methods 
have previously been suggested 0, as has detection of 
the spatially-dependent phase in an ionization scheme Q . 
The interference of condensates released from a double- 
well trap was the first evidence of their phase coher- 
ence. At a vortex core, the phase will always have a 
singularity which is independent of atomic and trap pa- 
rameters. Since there is coherence over most of the con- 
densates, we suggest that the phase slip occuring at a 
vortex core be used to detect it. When colliding a con- 
densate in the ground state with a (presumed) vortex 
state, the vorticity is measured by counting the number 
of fringes which are "skipped" as they cross the core. 

We briefly review the theory of the BEG far below the 
critical temperature beginning from the Gross-Pitaevski 
(GP) equation for the condensate wavefunction ip{x,t), 



i^--2V2v + F(x)V + -|VlV- 



(1) 



In the case of axisymmetric harmonic traps, the external 
potential is 



y(r,z) = i(r2 + AV) 



(2) 



Here the time is scaled in units of the reciprocal trap 
frequency, and all lengths are scaled in units of the 
transverse trap width. 



so that the dimensionless interaction strength is 



(3) 
(4) 

(5) 



where TV is the number of condensate atoms, and a is 
the s-wave scattering length. Since we are interested pri- 
marily in the two-dimensional dynamics, we assume that 
the z dimension is thin enough that all dynamics in that 
direction can be neglected (A » 1). 

We distinguish between pure vortex states and generic 
states with vorticity. The latter are any states for which 
there exists a closed curve C such that the circulation (in 
unsealed variables) obeys 
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Owing to the single-valuedness of the wavefunction, m 
must be an integer. The structure near the vortex core is 
always similar: for cylindrical coordinates (p, (f)) centered 
at (xc, yc)^ the density goes from to no within a radius 
of the order of the healing length 



(7) 



and the phase is 

(j) = ruLp = m arctan 



y-Vc 



(8) 



Girculation is conserved around any closed curve flowing 
with the condensate (Kelvin's theorem). On the other 
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hand, these states are not stationary; our numerical so- 
lution of the GP equation shows precession around the 
trap center with a period depending on the distance from 
it to the core and on u. For example, when u — 300, a 
core 1.5 from the center precesses with a period « 19 
(Fig. |l|). This appears to violate Ehrenfest's theorem; 
however, the precession is accompanied by considerable 
density excitations so that the center of mass still obeys 
classical equations of motion. From this we may expect 
that the dynamics of generic vorticity states depends not 
only on the atomic collisions but also damping rates of 
excitations. 

In the case of axisymmetric traps, a subset of the vor- 
ticity states are pure vortex states of the form 

V'(r,z,^)=/(r,z)e-'^ (9) 

These states are stationary and are eigenstates of az- 
imuthal angular momentum L^, = Hm. The mctastabil- 
ity of these states is still controversial [|l^,|ll|]. However, 
one decay scenario is that pure vortex states turn into 
vorticity states with the core off-center [|l^; the possible 
observation of this process is one of the motivations for 
our detection scheme. 

For interference between a generic vorticity state and 
a ground state moving towards each other with relative 
velocity 2k in the x-direction, in the core region the wave- 
function is approximated by 

V'(x) ~ VTITe'''^ + V^e(~''^^+™'^) (10) 

resulting in a density 

n(x) ~ rii + 712 + 2^/nin2 cos{2kx — rrup) (11) 

This density displays a series of fringes, with |m| -|- 1 
fringes merging into 1 at the origin, which is the location 
of the vortex core. Such patterns have been used to de- 
tect vortices in nonlinear optics fl^ . The momentum k 
may be acquired from free expansion of the condensates 
or by moving the traps. 

We have computed several interference patterns be- 
tween states with and without vorticity, by numerically 
solving the GP equation in two dimensions. Our initial 
configuration consists of a pair of trapped condensate 
clouds with equal numbers of atoms. It is not neces- 
sary that either the actual number or the relative phase 
between the clouds be known. In Fig. ||, the cloud cen- 
tered at (6.25, 0) begins in the m = — 1 pure vortex state 
while that at (—6.25, 0) is the trap ground state. The 
traps are switched off at t = 0. Up until t ~ 1.5 the 
two clouds expand independently. Then, as in the first 
interference experiment ||], the "pulsed" sources lead to 
straight fringes. However, as the core region of the upper 
condensate approaches the lower condensate, the differ- 
ence in velocity between the upper and lower sides of the 
vortex tilts the fringes from the vertical. Finally, as the 
fringes reach the position of the core, the phase slip in 
fringes appears clearly (Fig. ||). We note that the core is 
smaller than the fringe spacing in this example. 



In the second configuration, shown in Fig. the cloud 
on the right is replaced by an m = — 1 generic vorticity 
state with the trap center placed at (—12.5,0) and the 
core at {xc — 12.5, j/c = 0.75). This state was produced 
by evolving the wavefunction in imaginary time, subject 
to the constraints of normalization and that the phase 
be given by Eq. (||). The cloud on the left is the ground 
state centered at (12.5, 0). Also, the condensates are now 
launched towards each other so the wavefunction is mul- 
tiplied by exp(iA:a;) for x < and by exp(— ifca;) for a; > 0, 
with fc = 4. Since the vortex is off-center, it begins to pre- 
cess before the two clouds collide. In this example, the 
period of this precession is much longer than the time 
until the interference pattern reaches the core. Thus, in 
Fig. H the position of the phase-slip above the x-axis in- 
dicates the initial position of the vortex. This suggests 
that precession of a core may be observed by repeatedly 
creating an off-center vortex, allowing it to evolve in the 
trap for successively longer times, and then interfering it 
with a ground state. 

We have also examined the case of vortex states with 
|m| > 1. Any non-azimuthally symmetric perturbation 
is expected to break these into a collection of m = ±1 
generic vortices. Fig. ^ is the pattern resulting from 
freely expanding clouds at (—6.25,0), m = 0, and 
(6.25,0), m = 2. In this case two fringes are skipped; 
however it is not possible to determine from the density 
alone that in fact we have two nearby m = 1 phase slips. 

In an actual experiment, the core region itself may be 
obscured for several reasons. First, if the local density is 
high enough, the core radius ^ will be less than the imag- 
ing wavelength and may not be resolved. Second, light 
traveling along the z-direction will always pass through 
regions of non-zero density if the core is not exactly par- 
allel to the laser beam. (This problem can be avoided by 
optically pumping atoms in thin sheets perpendicular to 
the imaging beam [^). Finally, the position of the core 
is subject to thermal and zero-point fluctuations of the 
order of its radius ||lj]. We note that detection of the 
fringes avoids all of these problems. The same method 
may also be used to measure non-zero winding numbers 
in toroidal traps. 

We have assumed that damping takes place on longer 
timescales than the dynamics seen in our simulations, 
which run up to a few trap periods. As mentioned above, 
to accurately follow the motion of generic vorticity states 
requires that we go beyond the GP equation to include 
thermal effects. We are continuing research in this direc- 
tion. 
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FIG. 6. Density at i = 10 after free expansion of two 
trapped condensates, with right cloud starting in the m = 2 
vortex state, left cloud in the trap ground state. Interaction 
strength u = 300. 



FIG. 1. Position of a precessing core in a generic vorticity 
state. Dashed curve, x; solid curve, y coordinate. The initial 
state was prepared by imaginary time evolution subject to 
constrained normalization and phase. Interaction strength 
u = 300. 



FIG. 2. Initial density of two-trap BEC for interference. 
Right condensate is in the m = — 1 vortex state, left conden- 
sate is in the trap ground state. Interaction strength u = 300. 



FIG. 3. Density at t = 10 after the traps of Fig. g have 
been switched off. The original vortex is clearly evident from 
the missing fringe near (8,0). 



FIG. 4. Initial density for two condensates moving towards 
each other. Left cloud is the ground state, right cloud has an 
m — —1 vortex core at (—12.5,0.75). Interaction strength is 
u = 300. 



FIG. 5. Density at t = 10 after the traps of Fig. ^ moving 
at fc = ±4 have been switched off. The position of the phase 
slip is still above the a;-axis, and has precessed less than 1/8 
of the way around the cloud. 
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